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1 Introduction 



The theory of Quantum Chromodynamics (QCD) is a central component of the Standard Model of 
elementary particle physics. It accounts for the strong nuclear force in terms of the exchange of spin- 
one entities known as gluons; these carry a quantum number known as "colour" and are exchanged 
between themselves and between the fermions referred to as quarks, which carry electric as well as 
colour charge. There is a law of nature that the overall colour of a stable elementary particle must 
be zero. Thus the quarks and gluons are never observed alone, hnt mesons are formed (primarily) 
from quark-antiquark pairs and baryons from three quarks or antiquarks. The qiiarks that provide 
the basic constitution and quantum numbers of a hadronic particle are known as valence quarks. 

When a high-energy electron (or positron) collides with a proton, as in the HERA collider, the 
simplest nuclear reaction that can occur is that a quark is ejected to give a high-energy jet of particles. 
This is Deep Inelastic Scattering (DIS) and is mediated by the exchange, between the electron and 
the quark, of a virtual photon, W or Z. Scattering processes can involve valence quarks, but can also 
occur by means of virtual quark-antiquark pairs that appear temporarily within the proton. Indeed, 
the proton may be considered as consisting of a continuum of different combinations of partons, that 
is to say (anti)quarks and gluons, and depending on how violently it is struck by the exchanged boson. 
Thus the parton density functions (PDFs) of the proton are functions both of the fraction x of the 
protons taken by a given parton and of the virtuality of the exchanged boson. 

2 The HERA Collider 

The HERA Collider has been operating at DESY, Hamburg, since the year 1992. Electrons or positrons 
are stored at 27.5 GeV and collide with a stored proton beam of energy 920 GeV. The experiments 
HI and ZEUS are located at the collision points and record the various ep interactions that occur. 
Over the years, the collider has achieved increasingly high luminosities which have enabled the different 
types of process to be studied in greater depth and detail. During 2001-2002 an upgrade to the collider 
and experiments was carried out; however, earlier analyses have continued and many of the results 
presented here come from these. Using silicon and drift chamber systems, each of the experiments 
is able to record tracks of charged particles emanating from the interaction point, supplemented by 
electromagnetic and hadronic calorimeters to measure the overall energy of jets and high energy single 
particles. The apparatus is asymmetric, with calorimetry in the proton (forward, "positive") direction 
designed to record large hadronic energy deposition, and calorimetry in the electron ("negative") 
direction aimed at measuring precisely the scattered electron. Both electron and positron beams were 
used. 

3 Tests of QCD at HERA 

Tests of QCD at HERA may broadly be considered as of the following kinds: 

• The PDF's of the proton measured at different values should be consistent with each other, 
as described by the DGLAP equations. 

• Cross sections for processes with different numbers of QCD vertices should be predictable. Here, 

in particular, we are looking at processes with different numbers of high transverse momentum 
(pt) jets, many of which correspond to radiated gluons, each of which represents an additional 
QCD vertex. We note that the simple lowest-order (LO) DIS process has no vertices involving 
gluons, and is purely electromagnetic or electroweak in nature once the PDF is known. 

• The value of the strong coupling constant a^, which acts at quark-gluon and gluon-gluon vertices, 
must be consistent when measured in different ways. The value of varies ("runs") with the 
momentum exchanged at the relevant vertex, in a way specified by QCD theory. 
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Figure 1: Inclusive dijet production in DIS from ZEUS, with NLO QCD fit. 

Various types of process are useful for these studies. Owing to the complexity of the QCD calcula- 
tions, the processes that can be calculated theoretically are those that involve perturbative expansions 
in terms of as, starting at LO. Even next-to-leading order calculations (NLO) are highly complex, 
however. Many such calculations have now been completed, and a few NNLO calculations, with spe- 
cial methods (resummation) employed to take account of multiple initial-state gluon radiation from 
the proton system. Since the running value of as falls with the momentum exchanged in the collision 
process, high values of exchanged momentum are preferred, so-called "hard" processes. "Soft" pro- 
cesses require a higher-order or non-perturbative approach which is not usually achievable from first 
principles. The following types of hard process have been studied particularly at HERA; more than 
one of these features may of course be present in a given reaction. 

• DIS processes, i.e. those involving high of the exchanged boson. 

• Processes involving high-pr jets. 

• The production of heavy-quark systems, i.e. charm and beauty quarks. 

• Hard diffractive processes in which the exchanged colour-neutral object may be treated parton- 
ically. 

4 Jets in DIS 

A long-standing analysis program at HERA has involved studying the PDF's of the proton by means of 
the ejection of one or more jets by the incident virtual boson. At lowest order (LO) only a single parton 
is involved, and at a fixed value the process does not display explicit QCD features. However the 
variation of the PDF's with is given by QCD-governed evolution equations (the DGLAP equations). 
These have been well verified in processes dominated by u and d type quarks. When inclusive jets are 
measured, the contributions from multi-jet final states are included, and here QCD-governed processes 
contribute. Both HI and ZEUS have published QCD-based analyses of the parton structure of the 
proton. P 

The two experiments have recently measured inclusive jet production to improved precision in 
high-Q^ neutral-current DIS. It is convenient to perform DIS jet measurements in the Breit frame, in 
which the incoming virtual boson and proton define the event axis. If the outgoing state contains the 
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Figure 2: Results of fit for by HI to inclusive jet data, as a function of jet transverse energy. 

proton remnant and just one jet, the latter then has no transverse momentum relative to this axis. 
High-i?T jets in the Breit Frame therefore indicate a > 2-jet process which is sensitive to QCD physics. 
Dijet measurements by ZEUS are shown in Fig. ^ for a range of exchanged boson virtualities.[2I The 
jet cross sections have been evaluated as functions of several kinematic and jet variables and are well 
described by the NLO QCD calculation DISENT.p] The aim has been to identify regions of phase 
space where uncertainties from higher-order corrections are small, so as to facilitate the evaluation of 
improved gluon densities in the proton in global fits to the proton PDFs. 

HI have used the NLO QCD program NLOJET++ to represent their data at the parton level. |1] 
It is performed in the MS scheme with five flavours, the effects of hadronisation being evaluated using 
the generators RAPGAP (based on LO matrix elements) and DJANGOH (using the Colour Dipole 
Model). Used together with Q^-varying PDFs, this gives a good description of the data and allows 
a value of Og to be extracted by varying as in the fits. The fitted as values (Fig. [21) show running, with 
the equivalent value at Mz showing good consistency as a function of transverse jet energy. A value of 
0.01197 lb 0.0016 lb 0.0047 is obtained, in agreement with other determinations and with competitive 
precision. 

HI and ZEUS have measured the inclusive production of jets in DIS in the forward direction. 
Using DISENT, HI find that their cross sections (Fig. ^ are not well described by LO or NLO 
calculations in QCD. [21 Models with a resolved photon contribution, and the Colour Dipole Model 
(CDM) implemented in ARIADNE, give an improvement, [H| as does the CASCADE model which 
incorporates the CCFM description of the proton using /cT^-unintegrated parton densities. [Ij However 
none of the models considered is satisfactory, and HI conclude that a higher-order description is 
required. 

ZEUS however, without going to the very lowest x values attained by HI, find that DISENT is 
satisfactory (Fig. Q in an update of a previous analysis. ^Uj They assign large values to the QCD scale 
uncertainties in the perturbative calculation, and these cover the discrepancies between the latter and 
the data. It is evident that there is a theoretical issue here which needs to be resolved, and as yet 
these measurements do not give clear evidence in favour of the CCFM model. 

At high values of Q^, cross sections have been measured for charged-current DIS using the newly 
developed polarised beams at HERA (Fig. [HJ. The differences between the two cross sections 
as a function of polarisation are very evident and in agreement with the Standard Model, using the 
Monte Carlos ARIADNE and MEPS to describe the hadronisation effects. The earlier ZEUS NLO 
PDF analysis has been updated to include these data, and the resulting picture is entirely consistent. 
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Figure 3: HI measurements of jet production cross sections at low x in DIS compared to various 
theoretical models (see text). 
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Figure 4: Preliminary ZEUS measurements of forward jet production cross sections at low x in DIS 
compared to the DISENT model. 



5 Jet structure and event shapes 

The particles in a jet are produced by the soft fragmentation of the initiating high-p^^ parton. The 
QCD processes involved here are highly non-perturbative and must be modelled phenomenologically, 
as is done in the leading-logarithm Monte Carlos PYTHIA and HERWIG. However, harder radiation 
of partons within the jet can deform its shape from a long pencil-like object (with high thrust) into 
a form with broader and asymmetric characteristics. These can be described by various event-shape 
variables, whose values are given by an a^-dependent term together with so-called power corrections 
to account for the soft effects and ensure that these are suitably merged on to the hard processes. 
Both ZEUS and HI have extended their measurements of a variety of event-shape variables. jl21 [T^ 
and have extracted values of ag together with a parameter ao which, according to power-correction 
theory, should be universal with a value of ~ 0.5 (Fig. IH}. Differential distributions in the event shape 
variables have also been studied. The theoretical expectations are reasonably well observed in most 
of the event shapes measured, but there is a lack of consistency in the extracted qq and as variables 
that may indicate a need for higher order terms. The two experiments also disagree over the detailed 
values extracted from their fits, apparently because different kinematic regions were used. This gives 
an indication of the systematic uncertainties currently involved with the power correction method. 

A further probe of the hard structure of QCD jets comes from the study of subjets. These 
are defined as jet-like substructures within jets, which may be obtained by reapplying the jet cluster 
algorithm at a smaller resolution scale ycut than was used for identifying the jet itself. Sub jet have been 
studied by ZEUS in neutral current DIS. The jets were identified using the kx cluster algorithm in the 
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Figure 5: (Left) charged-current scattering cross sections of electrons (black) and positrons (grey) 
from ZEUS and HI as a function of lepton polarization, compared to standard model predictions. 
(Right) Parton densities from the ZEUS fits to jet data (ZEUS- JETS), and incorporating the new 
polarisation data (ZEUS-poI); the uncertainties are those on the ZEUS-pol fits. 
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Figure 6: Measurements by HI and ZEUS of event shape variables: r, T = thrust; c-parameter; B = 
broadening; M'^,pq =invariant jet mass. Subscripts T, C denote the thrust axis (ZEUS, HI), 7 the 
incoming photon axis (ZEUS); unsubscripted r and B (HI) are relative to the photon axis. 



longitudinally invariant inclusive mode in the laboratory frame. Measurements of subjet distributions 
were carried out as functions of various kinematic quantities, including the ratio of the subjet transverse 
energy to that of the jet, and the differences between the pseudorapidities and the azimuth values of 
the subjet and of the jet. The measured distributions are used to study the pattern of parton radiation 
by comparing them with Monte Carlo models and perturbative QCD calculations. It is found that 
ARIADNE performs better than LEPTO-MEPS and that comparison with NLO predictions allows a 
confirmation that quark rather than gluon jets are mostly being observed (Fig. [71). 



6 Multi-jet properties 

When two or more jets emerge in a high energy process, a QCD-based description is essential even 
at lowest order. HI have measured the correlations in the azimuthal angle when measuring two-jet 
final states in DIS. This distribution is sensitive to additional parton radiation both in the initial state 
and from the outgoing hard partons. HI find that the process is sensitive to higher order radiative 
processes and the simple NLO two-jet amplitudes are insufficient to describe the data. NLOJET three- 
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Figure 7: Kinematic properties of subjets in DIS from ZEUS: energy fraction (left) and pseudorapidity 
difference (right), compared to NLO prediction using DISENT. 
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Figure 8: Azimuthal correlations between jets in two-jet final states in DIS measured in the hadronic 
centre-of-mass frame. The NLO predictions are normalised to the data outside the highest A0* bin, 
while the CASCADE predictions are absolute. 



jet calculations give improvement but are still insufficient (Fig. IHta)), as are attempts to incorporate 
a resolved-photon component in the RAPGAP model. The CASCADE model also fails, suggesting 
a need for more theoretical input. In principle, data of this kind can give information on the kT- 
unintegrated parton densities within the proton. 

Similar measurements from ZEUS, however, indicate good agreement with NLO QCD calculations 
for the dependence of the overall cross sections on kinematic variables such asQ'^, x and jet .E'r.|14j Like 
HI, ZEUS require higher order terms in order to describe the azimuthal correlations. The inclusion of 
O(a^) contributions increases the predictions by an order of magnitude when the two highest- jets 
are not opposite in azimuth. Overall, the description of the two- and three-jet processes is reasonable 
although not perfect at the highest jet-E'r values. 

7 Prompt photon production 

Prompt photons are those that emerge directly from a hard scattering process, as opposed to photons 
found as products of decay processes or as soft radiation from an outgoing particle. They can give 
information about QCD processes in a more direct way than do jets, which require hadronisation 
effects to be allowed for in order to be interpreted as a measure of quark and gluon properties. 
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Figure 9: Distributions of prompt photons in DIS from HI, as a function of photon transverse energy 
and pseudorapidity. The theoretical curves LL and QQ show the contributions of radiation off the 
electron and the quark line respectively, the interference term being small. 
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Figure 10: ZEUS prompt photon measurements in conjunction with a jet: a) differential photon cross 
sections as a function of transverse energy, b) as a function of laboratory pseudorapidity. 

Following a first publication by ZEUS on prompt photons in DIS, HI have presented differential 
DIS cross sections. jl5j They are obtained for photon transverse energies of 3 - 10 GeV and labora- 
tory pseudorapidities of —1.2 - 1.8, and the data are compared with the predictions of a new LO 
calculation 16 , which gives reasonable agreement. Comparisons with the predictions of the event 
generators PYTHIA and HERWIG are also presented. Fig. IHl illustrates these results. 

In photoproduction, the production of prompt photons with an accompanying jet has been studied 
by ZEUS, using a different technique from their earlier results. [T7| In the present case the photons 
were identified from their conversion probabilities in a preshower detector installed in front of the 
calorimeter. The differential 7 + jet cross sections were reconstructed as functions of the transverse 
energy, pseudorapidity and x-y, the fraction of the incoming photon momentum taken by the photon- 
jet system. Predictions based on leading-logarithm parton-shower Monte Carlo models were compared 
with the data (Fig. llUj) and were found generally to underestimate the cross sections for photons of 
transverse energy below 7 GeV, while the NLO QCD calculations (applying the jet algorithm to the 
partons) agree with the data better. For higher transverse energies of the prompt photons both types 
of calculation are in good agreement with the data. 

8 Jets in photoproduction 

HI have measured the photoproduction of dijet systems with improved statistics fFig. Illjl.'lS' When 
the dijet mass is above 65 GeV, the distributions in | cos 9*\ show a clear difference between the direct- 
and resolved-dominated samples. In direct processes, the incoming photon interacts with the proton 
by coupling to a high-ET quark line: the hard exchanged object is thus a virtual quark, with spin 
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Figure 11: Dijet production in photoproduction by HI: cross sections for resolved-dominated events 
(left) and direct-dominated events (right) as a function of dijet centre-of-mass scattering angle. 
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Figure 12: Three and four jet final states in photoproduction in ZEUS: the fraction of photon mo- 
mentum carried by the jets, compared to some theoretical models. 



one-half. In resolved processes, a parton from the partonic substructure of the photon interacts with 
a quark or a gluon from the proton, and the exchanged object is more likely to be a spin-one gluon. 
Hence the QCD dynamics is different: a less steep | cos0*| distribution occurs in the direct case owing 
to the exclusively spin-half propagator. 

The data are well described by PYTHIA, in which the QCD dynamics is of course inbuilt, rescaled 
by a factor of 1.2; this does not affect the shape of the distributions. An NLO calculation describes a 
variety of distributions to within in%.|19j 

Three- and four-jet final states have been measured in photoproduction at HERA by ZEUS. The 
integrated luminosity of 121 pb~^ represents over seven times the luminosity of the previous HERA 
publication on three jets in photoproduction, while the four-jet photoproduction cross section has been 
measured for the first time and constitutes the highest order process so far studied at HERA. The 
events have been studied for masses of the jet system greater than 25 GeV and greater than 50 GeV. 
A comparison with the PYTHIA and HERWIG shows that a multiple parton interaction scheme is 
necessary, but that tuning to ep data is necessary; this has been done with HERWIG but not yet with 
PYTHIA (Fig. 1121) . The three-jet cross sections have been compared to an O(aa^) perturbative QCD 
calculation. |2nj This describes the higher-mass region reasonably well but underestimates the data for 
jet system masses below 50 GeV. 
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Figure 13: The fraction of the photon momentum taken by a dijet system containing charm or beauty. 
Events are selected with two or more jets of transverse momentum above 11, 8 GeV respectively in 
the central laboratory pseudo-rapidity range —0.9 - 1.3. 



9 Heavy flavour studies 

The dominant process for heavy flavour production at HERA is photon-gluon fusion: jg — > QQ. The 
cross sections therefore depend strongly on the gluon density in the proton. There are two main 
theoretical approaches, namely the zero-mass variable flavour number scheme (ZM-VFNS) which is 
suitable for values much higher than the quark mass squared, and the fixed-flavour number scheme 
which allows for quark masses, but which in present implementations becomes inaccurate at high 
owing to large terms of the form logQ^/M^ in the perturbative expansion. Combinations of the two 
schemes have been devised, referred to simply as VFNS. 

An analysis of charm and beauty dijet photoproduction cross sections has been carried out by 
Hl.|21j The fractions of events containing charm and beauty quarks are determined using a method 
based on the impact parameter in the transverse plane of tracks relative to the primary vertex, as 
measured by the HI central vertex detector. A variety of differential dijet cross sections for charm 
and beauty and their relative contributions to the flavour inclusive dijet have been calculated, as 
illustrated in Fig. 113(1 . Comparison is made with an NLO QCD calculation. j23j Taking into account 
the theoretical uncertainties, the charm cross sections are consistent with a QCD calculation at NLO, 
while the predicted cross sections for beauty production are somewhat lower than the measured values 
in the resolved-dominated region of the distribution. For direct-dominated events, the respective 
contributions of charm and beauty events to the total dijet cross section should be 4/11 and 1/11, 
given by the squares of the charges of the respective quarks, for dijet masses much greater than the 
quark mass. This is found to hold within experimental uncertainties. 

These results follow a measurement of the beauty production cross section in ep collisions, both 
photoproduction and DIS, in which events were selected by requiring the presence of jets and muons 
in the final state. ^2] Both the long lifetime and the large mass of b-flavoured hadrons were exploited 
to identify events containing beauty quarks. The differential cross sections in photoproduction, and 
in deep inelastic scattering, with 2 < < 100 GeV^ were compared with perturbative LO and NLO 
QCD calculations. As evidenced in Fig. 1141 the predictions are found to be somewhat lower than the 
data. 

Beauty photoproduction with events in which two muons are observed in the final state has been 
measured with the ZEUS detector at HERA (Fig. I14| right). A low pT threshold for muon identi- 
fication, in combination with the large rapidity coverage of the ZEUS muon system, gives access to 
essentially the full phase space for beauty production. The dimuon selection suppresses backgrounds 
from charm and light flavour production. The total cross section for beauty production has been 
measured and differential cross sections for b quark production have been extracted on the basis of 
our knowledge of b fragmentation. A comparison to NLO QCD predictions j23j and to other ZEUS 
measurements is shown in Fig. 1141 Although no single measurement is significantly different from 
theory, taking theoretical and experimental uncertainties into account, there is a consistent trend for 
theory to be low at lower jet momenta, although high and low values of do not appear to differ in 
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Figure 14: (Left) HI cross section for b quark production in DIS as a function of compared to 
theoretical models. (Right) cross section for b quark photoproduction in ZEUS (preliminary), with 
earlier results and NLO QCD prediction. 



this respect. 

HI have presented new results for the proton structure function measured for the production of 
final states involving c and b quarks (Fig. I15|) _24_. Beauty and charm comprise respectively 0.4 - 3% 
and ~24% of the total cross section over this kinematic range. These are the first such measurements 
for b quarks and were made exploiting the HI silicon vertex detector. The slopes of the distributions 
with at given x (scaling variations) arise from perturbative-QCD effects on the PDF's and are in 
accord with predictions. 

The HI data are compared with several VFNS-based theoretical approaches incorporating different 
proton PDFs. The agreement is reasonable for both charm and beauty, but the present experimental 
and theoretical uncertainties are such that the different models cannot yet be distinguished. 



10 Diffract ive production 

Some recent results on diffractive production in DIS are illustrated in Fig. 1161 The quantity cjt?^^^ 
measures the inclusive differential cross section and is a good approximation to the diffractive proton 
structure function. The quantities xjp and f3 approximate respectively to the fraction of the pro- 
ton momentum carried by the colourless exchanged "pomeron", and to the fraction of the latter's 
momentum carried by the struck quark. 

The figure shows results obtained by HI at xjp = 0.002 for two values of Q^. There is good 
agreement with the predictions of a model in which the "pomeron" is modelled as two gluons.[25] The 
results are also in good agreement with earlier HI and ZEUS diffractive measurements. j^Hl 



11 Conclusions 

The results presented here have of necessity been selective, chosen to illustrate the range and depth of 
QCD studies that are being carried out at HERA. Much progress has been made over recent years, in 
the type of studies that can be performed, the precision achieved, and in theoretical understanding. 
HERA will complete its data-taking phase in mid-2007. The final analysis of its results will provide a 
unique corpus of knowledge that will not be overtaken for the foreseeable future. 
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Figure 15: Measurements of the proton structure function F2 by HI for processes involving production 
of h and c quarks in the final state, compared to models with different proton PDFs. 
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